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Abstract A comprehensive understanding of the solar magnetic cycle requires de- 
tailed modeling of the solar interior including the maintenance and variation of large 
scale flows (differential rotation and meridional flow), the solar dynamo and the flux 
emergence process connecting the magnetic field in the solar convection zone with 
magnetic field in the photosphere and above. Due to the vast range of time and 
length scales encountered, a single model of the entire convection zone is still out of 
reach. However, a variety of aspects can be modeled through a combined approach 
of 3D MHD models and simplified descriptions. We will briefly review our current 
theoretical understanding of these processes based on numerical models of the solar 
interior. 



1 Introduction 

The solar convection zone comprises the outer most 30 % of the solar radius and 
contains about 2 % of the total solar mass. Due to a density variation of more than 6 
orders of magnitude a variety of different physical regimes are encountered. While 
fluid motions are are highly subsonic (Ma « 10~ 4 ) and strongly influenced by rota- 
tion near the base of the convection zone, they turn supersonic in the photosphere 
and the influence of rotation diminishes. The pressure scale height varies between 
about 50 Mm at the base of the convection one and about 100 km in the photosphere 
of the sun. As a consequence a comprehensive model of the entire convection zone 
is currently out of reach and different aspects have to be modeled independently. 
The deep convection zone up to about 10-20 Mm beneath the photosphere can be 
modeled most efficiently using the anelastic approach which is filteri ng out sound 



waves, but is fully considering the compressibility in the stratification (iGlatzmaiei 
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19841) . T he upper most parts of the convection zone require fully compressible MHD 



(see e.g. lNordlund et al.L 120091 for a recent review). While most anelastic models of 



the solar interior are global models with computational domains covering an en- 
tire shell between two radii (or at least a shell segment), MHD models of the outer 
parts of the convection zone typically focus on details in rectangular computational 
domains. 

Apart from 3D MHD models adapted to the different physical regimes a vari- 
ety aspects have been modeled based on simplified models, such as the mean field 
approach. Here the focus is on the large scales, while the effects of unresolved tur- 
bulence is parametrized. Non-linear terms in the momentum, energy and induction 
equations lead to non-vanishing second order correlation terms of small scale quan- 
tities that act as drivers for large scale flows or as turbulent induction effects for the 
large scale magnetic field. The decomposition into large and small scale properties 
and the arising correlation terms driving large scale flows are the strength and the 
weakness of this approach at the same time. On one hand the computational ex- 
pense is decreased by orders of magnitude allowing for simulations covering long 
time scales as well as exploring wide parameter ranges, on the other hand the results 
are heavily dependent on parametrization of the second orde r correlation terms. For 
a comp rehensive description of mean field theory we refer to Riidiger & Hollerbachl 
d2004. 



2 Differential rotation and meridional flow 



Differential rotation is the consequence of angular momentum transport in the solar 
convection zone. Starting with a decomposition of the turbulent velocity field into 
fluctuating and (axisymmetric) mean flows v = (v) + V leads to the following terms 
in the angular momentum flux (neglecting magnetic terms for simplicity): 



(/v}=rsin0p( (v/v/) + (v,)i2rsine ) 
Reynolds stress Meridional flow 



(1) 



Here the Reynolds stress describes the transport due to correlations of fluctuat- 
ing turbulent velocity components, while the Meridional flow describes the trans- 
port due to large scale coherent mean flows in the r — 6 plane. Angular mo- 
mentum transport through Reynolds stresses requires the presence of rotation and 
anisotropy and expressi ons for these transport te r ms have be e n deri ved within the 
mean fiel d approach by iDurnev & Spruitl ( 1979b : Hathaway ( 1984 ) and more re- 
cently by iKitchatinov & Riidiger! ( 1993b using a quasi-linear approach (see also 
Kitchatinov & Riidiger d2005h for an improved representation). In 3D simulations 
the influence of rotation on convection leads to a preferential north-south align- 
ment of convection cells ( Gilman , 1979 ; Miesch et al.ll2000llBrun & Toomrell2002l : 



Miesch et al.L |2008). The consequence is a dominance of east-west motions over 



north-south motions. By means of the Coriolis force eastward (faster rotating) flows 
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are deflected equatorward, while westward (slower rotating) flows are deflected 
poleward, leading on average to an equatorward transport of angular momentum. 




Fig. 1 Contour plots of differential rotation (a), entropy perturb ation (b) and stream function of 
meridional flow (c) using the mean field model o f lRempelll2005h . Panel (d) shows the differential 
rotation profile obtained using the same parametrization of the Reynolds stress but neglecting the 
effects of baroclinicity. 



The profile of differential rotation cannot be determined on the basis of angular 
momentum transport processes alone. As stationary state requires beside vanishing 
divergence of the total angular momentum flux also a force balance in the meridional 
plane between Coriolis, centrifugal, buoyancy and pressure forces. The latter is most 
conveniently expressed by (follows from ^-component of vorticity equation): 

dQ. 2 g ds 
rsm9— = — (2) 
dz c„rd& 
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Helioseismic observations bv lThompson et al. (1996) show clearly a differential ro- 
tation profile with contours of constant £2 inclined by about 25 deg with respect to 
the rotation axis (deviation from Taylor-Proudman state). It turns out that avoid- 
ing the Taylor-Proudman state is a key problem for a theoretical understanding 
of solar differential rotation. While early models atte mpted to achieve th i s by a s- 



suming large viscositie s ('Tay lor-number puzzle' after iBrandenburg et al.1 ( 19901) ). 



iKitchatinov & Riidigei (1995) showed that an alternative solution of this problem 
can be given if the anisotropic convective energy transport is considered, leading to 
a pole-equator temperature difference of about 10 K. Anisotropic convective energy 
transport is automatically considered in global 3D simulations, but in many cases it 
turns out to be insuff i cient f or obtaining solar-like differential rotation. 

Recently Rempell (2005) showed that coupling between the tachocline and con- 
vection zone can also provide the latitudinal entropy variation needed to explain the 
observed profile of solar differential rotation. A typical solution from that model is 
shown in Fig. <£TJ, displaying differential rotation (a), corresponding entropy pertur- 
bation (b) and the stream function of the meridional flow (c). Panel (d) shows for 
comparison the profile of differential rotation obtained if the effects of the entropy 
perturbation displayed in (b) are neglected. Inclusion of this effect through the bot- 
tom boundary condition in 3D m odels allows also fo r more solar-like differential 



rotation in 3D convection models dMiesch et al. ■ 120061) 



While there is a general agreement that thermal effects are essential for solar- 
like differential rotation, it is still unclear whether the required latitudinal entropy 
variation is a consequence of anisotropic convective energy transport, imposed by 
the tachocline, or a combination of both. 

How does the meridional flow come into play here? A stationary solution re- 
quires that the divergence of the angular momentum flux Eq. (Q]) vanishes. While 
in very special situations the Reynolds stress can be divergence free on its own, in 
general a contribution from the meridional flow is required to close the system. It 
turns out that primarily the component of the Reynolds stress that transports angu- 
lar momentum parallel to the axis of rotation influences most strongly the direction 
of the meridional flow. If the transport of angular momentum is inward directed, 
the resulting meridional flow is poleward at the surface and equatorward near the 
bottom of the convection zone. While this is found in most mean field models such 



as IKitchatinov & Riidigerl (11995b . 3D simul ations present a more c omplicated sit- 
uation. Early models with lower resolution (iBrun & Toomrel 120021) typically show 
multi-cellular flow pattern, while a recent high resolution run dMiesch et al. , 2008 ) 
shows a single flow cell (poleward at top, equatorward at bottom of CZ) in the bulk 
of the convection zone. To which degree these results are converged with respect to 
numerical resolutions remains to be seen in the future. 

Differential rotation shows also cyclic variations known as torsional oscillations, 
which point towar d a close relation to the solar magnetic cycle. We refer here to 



Howel (120091) and IBrun & Rempell (120091) for reviews of observations as well as 



theoretical models for the time varying component of solar differential rotation. 
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Similar to models of differential rotation and meridional flow we discussed in sec- 
tion^ also the solar dynamo is modeled through a combination of mean field mod- 
els and 3D simulations. Currently mean field models of the solar dynamo are the 
only models that provide dynamo solutions that are compatible with basic cycle 
features and can be evolved over time scales much longer than a cycle. However, as 
already stated above, these models are heavily dependent on parametrization of tur- 
bulent induction processes and cannot provide an explanation from first principles. 
On the other hand, 3D MHD simulations describe currently only aspects of the dy- 
namo process, a comprehensive model of a solar dynamo with features compatible 
with the basic dynamo constraints derived from the solar butterfly diagram is still 
an open challenge. 

Regardless of the adopted modeling approach the primary uncertainties regard- 
ing the underlying dynamo process are similar. Many of these uncertainties result di- 
rectly from our limited ability to model processes from first principles and the rather 
sparse observational constraints on the solar interior. The best known ingredient is 
differential rotation (Q -effect) due to observational constra ints from hel ioseismol- 
ogy on the mean profile and variation of differential rotation ( IHowe , l2009h . But even 
the exact knowledge of the differential rotation profile is not sufficient to determine 
whether radial velocity gradients at the base of the convection zone (tachocline) 
or latitudinal gradients in the bulk of the convection zone play the major role in the 
generation of toroidal magnetic fields, since this would require knowledge of the de- 
tailed distribution of poloidal field in the convection zone. Even less known are the 
processes related to the regeneration of poloidal field (a-effect). A third unknown 
are the transport processes of magnetic flux in the convection zone. Since in general 
the locations where the a-effect and Q -effect operate do not coincide, transport of 
magnetic flux inbetween these regions is crucial for a coherent operation of the large 
dynamo in the convection zone. 

It is beyond the scope of this paper to review all the possib l e dyn amo scenar- 
ios which have been considered and we refer to ICharbonneau ( 2005 ) for further 
reading. In the following three subsections we will briefly discuss some of the key 
uncertainties and open questions. 



3.1 Role of tachocline 



Soon after helioseis mology revealed t he de tailed structure of differential rotation 



in the solar interior ( [Thompson et aU 119961) it was suggested that the base of the 



convection zone with its strong radial shear layer (tachocline) is a likely location 
for the solar dynamo (production of strong toroidal field by shear). In addition, the 
stable stratification found in the solar overshoot region at the base of the convection 
zone allows for storage of magnetic field over time scales comparable to the solar 
cycle. Both aspects are crucial since simulations of rising flux in the convection zone 
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as well as studies of magnetic stability in the solar overshoot regions (see section|4] 
for further detail) point toward a rather strong toroidal magnetic field of 10 5 Gauss 
at the base of the convection zone. More recently the role of the tachocline for the 
organization and amplification of large scale tor oidal field has been a lso seen in 
global 3D MHD sim ulations of the solar dynamo ([Browning et al. ■ 120061) . However, 



Brown et a l. (2009) presented simulations of solar like stars at faster rotation rates, 



which point toward the possibility that substantial magnetic field can be produced 
and maintained within the convection zone in near equatorial regions. It is currently 
not clear to which degree this result can be also relevant for the solar rotation rate. 

While most models of flux emergence point toward a field strength of 10 Gauss 
at the base of the convection zone, it is far from trivial to amplify field to this strength 
solely through diff erential rotation . Dynamo models that include non-linear feed- 
backs consistently (IRempell 120061) lead to an upper limit more ar ound 10 4 Gauss, 



simila r values are also found in most 3D simulations such as dBrowning et al 



20061) . Whether this discrepancy can be bridged through an alternative field amplifi- 



cat ion mechanism (e.g. harvest ing potential energy of the stratification as proposed 
by iRempel & Schiisslerl (120011) ^ or the possibility that also initially weaker mag- 
netic field from the bulk of the convection zone can lead to active region formation 
is currently an open question. 



3.2 Regeneration process ofpoloidal field 



The details of the processes rebuilding the poloidal field from toroidal field are still 
very uncertain. In the meanfield language these processes are formally described 
as a-effect and in the context of solar dynamo models the following 3 classes of 
a-effects are typically considered: 1. Helical turbulence, 2. MHD shear flow in- 
stabilities in the tachocline, 3. Rising flux tubes in the convection zone (Babcock- 
Leighton). While all these processes are likely to contribute, their amplitude and 
spacial distribution is not known well enough to clearly quantify their individual 
role. 

Furthermore recent research also points toward highly non-linear and also time 
dependent a-effec ts resulting from additional c onstraints due to conservation of 
magnetic helicity ( Brandenburg & Sandirl 2004 ). Indirect constraints on the oper- 
ation of the a-effect might be gained from helicity fluxes observable in the photo- 
sphere and above. 

The only a-effect contribution that is directly constrainable thro ugh observations 
is the Babcock-Leighton a- effect dBabcockl ll96lLlLei ghtqnl 1 1 969h. which has been 
used i n most of the recent flux-transport dynamo models dDikpati & Charbonneau , 
19991: iDikpati et all 12004 IRempell 120061) . The Babcock-Leighton a-effect is based 
on the flux emergence process leading to the formation of active regions, the key in- 
gredient is the systematic tilt resulting from the action of the Coriolis force twisting 
the rising flux tube. While it is possible to construct dynamo models entirely based 
on the Babcock-Leighton a-effect, these models lead in general to rather strong po- 
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lar fields at the surface in contradiction with observations, unless a strong magnetic 
diffusivity gradient and additio nal contributions fro m a-effects at the base of the 



convection zone are considered (IDikpati et all 120041) 



3.3 Transport of magnetic flux in convection zone 



Traditionally most models considered only turbulent transport in the convection 
zone, which can be decomposed (in the meanfield language) into diffusive transport 
(turbulent diffusion) but also advection like transport in form of turbulent pump- 
ing. The latter has been also studied extensively through 3D MHD simulations 



dTobias et al.L 11998. 2001). If magnetic field becomes sufficiently strong magnetic 



buoyancy drives additional transport in terms of rising flux bundles that can lead to 
the formation of active regions on the visible surface (see section|4]for more detail). 
Additional to these processes magnetic flux can be transported by the large scale 
meridional flow in the convection zone. Dynamo models based pri marily on the lat 



ter are call ed flux transport dynamos and were first introduced by IChoudhuri et al 



( 1995 ) and Dikpati & Charbonneau ( 19991) and have been developed further by sev- 
eral groups since then. 

The attraction of flux transport dynamos comes primarily from the fact that a 
meridional flow which is poleward at the top and equatorward at the bottom of the 
convection zone gives a very robust explanation for the equatorward propagation 
of the activity in the course of the solar cycle. In addition the poleward flow in the 
near surface levels in combination with the systematic tilt angle of sunspot groups 
leads automatically to the c orrect phase relati on between toroidal and poloidal field. 
However, as pointed out bv lSchusslerl ( 12005b . the phase relation is primarily a con- 
sequence of the tilt angle of active regions and in that sense only a weak constraint 
on dynamo processes in the solar interior. For meridional flow speeds consistent 
with surface observations and an extrapolation based on mass conservation for the 
deeper layers, these models also yield dynamo periods in agreement with the solar 
cycle. 

Overall the flux transport picture is currently one of the most successful sce- 
narios for the large scale solar dynamo, but (as many other models) it is based on 
two strong assumptions which cannot be proven from first principles: 1 . The merid- 
ional flow is dominated by one flow cell with poleward flow close to the surface 
layers and equatorward flow at the base of the convection zone. 2. Turbulent trans- 
port processes are sufficiently weak to allow advection effects to dominate. While 
meanfield models of differenti al rotation and meridion al flow typically lead to the 
required flow patterns (see e.g. lKiiker & Riidi gen 120051) . the situation is more com- 
plicated in 3D simulati ons. Most of the earlier models at moderate resolution lead 
to multi-cellular flows (IBrun & Toomrei 120021) . in contrast a more recent model at 
high er resolution shows a flow pattern dominated by a single cell in the convec- 



tion dMiesch et all 120081) . Overall the situation cannot be considered as converged 



yet. The amplitude of turbulent transport estimated from simple mixing length ar- 
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guments is typically 1-2 orders of magnitude larger than the values required for the 
flux transport picture. Since turbulent transport is in general more complicated than 
a simple diffusive transport this aspect needs to be studied in more detail through 
3D simulations taking into account the presence of large scale flows and the full 
non-linearity of the problem . 



4 Flux emergence process 



It is generally accepted that sunspots form from magnet ic field rising from the 
base of the convection zone to the surface (see reviews by iMoreno-Insertisi 1 19971 : 
Fisher et al. ■ l2000l:lFanll2004l. and further references therein), Solar dynamo mod- 
els as presented in section[3]focus on the large scale evolution of magnetic field an 
cannot address detailed processes such as the flux emergence process leading to the 
formation of sunspots on the visible surface of the sun. The latter is primarily a con- 
sequence of limited numerical resolution. Nevertheless, studying flux emergence is 
integral to our understanding of solar magnetism, since it allows us to connect the 
solar dynamo to observational constraints on the magnetic field structure in the solar 
photosphere. To date the flux emergence process has been studied decoupled from 
dynamo models using a variety of different approaches. 



4.1 Flux emergence in lower convection zone 



Early work was based on the thin flux tube approximation (IChoudh uri & Gilman, 



19871: 



Fan et al 



1993L 1 1994 iMoreno-Insertis et all 1 19941 ISchiissler et all 1 1994 



Caligari et all Il995l) . These studies concluded that the overall properties of active 



regions, such as the low latitude of emergence, latitudinal trend in tilt angles as well 
as asymmetries between leading and following spots can be understood on the basis 
of rising flux tubes, provided the initial field strength at the base of the convection 
zone is around 100 kG. This conclusion was also c onsistent with stability consid 



erations of flux in a subadiabatic overshoot region (iFerriz-Mas & Schiisslerl. 11993 
1995). 



Bas ed on two-dimensional MHD simulations it was early realized bv lSchiisslei 
19791) that untwist ed magnetic flux tubes cannot rise coherently and fragment. It 



was sh own later by IMoreno-Insertis & Emonet dl996h : lEmonet & Moreno-Insertis 



(1199a) that this fragmentation can be alleviated provided that flux tubes have enough 
initial twist. 

More recently also 3D MHD simulation s of rising fl ux tubes based on the anelas- 
tic approximation have become possible (IFanl [20083 and give support for results 
from earl ier simulati ons based on the thin flux tube approximation. It was however 
found by IFanl (12008b that there is a very delicate balance between the amount of 
twist required for a coherent rise and the amount of twist allowed to be in agree- 
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ment with observations of sunspot tilt angles (twist with the observed sign produces 
a tilt opposite to the effect of Coriolis forces on rising tubes). 

The simulations presented above consider the flux emergence process decoupled 
from convection. First attempts to address flu x emergence in globa l simul ations of 



the convection zone were made recently by iJouve & Brunl d2007 , l2009h . Under 



standing the interaction of emerging flux with the ambient convective motions in 
the convective envelope is a crucial step toward more realism; however, currently 
the focus on the global scale limits the resolution required to resolve this interaction 
in detail. Substantial progress will likely happen in the next decade with increase in 
computing power. 



4.2 Flux emergence in upper convection zone 

Another challenge is understanding the last stages of the flux emergence process in 
the near surface layers (last 10 - 20 Mm). All of the models presented above exclude 
the uppermost 10-20 Mm since the adopted approximations (thin flux tube, anelas- 
tic) loose their validity and also the steep decrease of pressure scale height and in- 
crease in convective velocities would lead to very stringent resolution and time step 
constraints. The upper most layers of the convection zone require fully compressible 
MHD and also the inclusion of radiative processes is necessary if a detailed compar- 



ison with th e available high resolution observations is desired (ICheung et all 12007 



2008 , 2010l) . While the primary modeling focus in the deep convection zone lies on 



large scale properties of active regions, near surface simulations focus on the de- 
tailed interaction of emerging flux with convection on the scale of granulation and 
below. One of the major open questions concerns the re-amplification process of 
magnetic flux into coherent sunspots from flux that has risen through a convection 
zone with a density stratification of six orders of magnitude. 

Recently MHD simulations with radiative transfer also provided a breakthrough 
in our understanding of sunspot fine structure in the photosphere such as umbral 
dots, penumbral filaments, light bridges and the Evershed flow in terms of a com- 
mon magneto-convection process modulated by inc l ination angle and field strength 
dSchussler & VoglerlEooalHeinemann et all 120071: iRempel et al.Ll2009blfah . 



4.3 Open questions, connection to dynamo models 

While we have seen strong progress in modeling the flux emergence process over 
the past decades, we do not have at this point a fully consistent model. The latter 
is a consequence of the fact that different aspects are modeled independently due 
to computational constraints. As a consequence there are some 'missing links' be- 
tween different models, which have to be addressed in the future through a more 
coherent coupling of models. Here we mention just a few of the open questions: 
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1 . As pointed out before most models of emerging magnetic flux require an initial 
field strength of about 10 s Gauss at the base of the convection zone to be consistent 
with observational constraints. On the other hand the majority of dynamo models 
falls short of such values, more typical are 10 4 Gauss. 2. Due to the strong density 
stratification in the convection zone even magnetic field with initially 10 5 is weak- 
ened to sub kG field strength in the upper most layers of the convection during the 
emergence process. It is currently an open question if such weak field can get re- 
amplified to sunspot field strength. Near surfac e simulations start ve ry often from 10 
kG field about 5 Mm beneath the photosphere dCheung et aUl2008l) to overcome the 
influence of convective motions. 3. Rising i2-shaped flux tubes in the deep convec- 
tion zone form typically as low wave number instability (m = 1 and m = 2 modes 
are preferred). In the near surface layers such low wave numbers should lead to 
much faster diverging motions in bipolar sunspot groups as observed (due to mag- 
netic tension) if sunspots stay connected to their deep roots. A possible dynam ical 
'disconnection' mechanism has been proposed by ISchiissler & Rempell ( 12005b . but 
it is also unclear if sunspots are sufficiently stable if they are rather shallow. 



5 Summary 

We presented here a brief summary of the state of the art of modeling of dynamical 
processes in the solar convection zone with focus on differential rotation/meridional 
flow, the large scale solar dynamo and the flux emergence process. We see currently 
in this field a dramatic change from more simplified models toward large 3D MHD 
simulations, primarily driven by the strong increase in computing resources. At the 
same time the field suffers from very limited observational constraints on processes 
in the solar interior. Progress in the future will heavily rely on improving and ex- 
ploiting helioseismic constraints and also on coupling models to allow for a check 
of consistency. In the near terms the latter is likely to be most successful for models 
of the flux emergence process. 
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